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Compact quantum electrodynamics in 2 + 1 dimensions often arises as an effective theory for a 
Mott insulator, with the Dirac fermions representing the low-energy spinous. An important and 
controversial issue in this context is whether a deconfinement transition takes place. We perform a 
renormalization group analysis to show that deconfinement occurs when N > Nc = ^ 1.161, 

where A'^ is the number of fermion replica. For N < Nc, however, there are two stable fixed points 
separated by a line containing a unstable non-trivial fixed point: a fixed point corresponding to 
the scaling limit of the non-compact theory, and another one governing the scaling behavior of the 
compact theory. The string tension associated to the confining interspinon potential is shown to 
exhibit a universal jump as A . Our results imply the stability of a spin liquid at the physical 

value N = 2 for Mott insulators. 
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An important topic currently under discussion in con- 
densed matter physics community is the emergence of 
deconfined quantum critical points in gauge theories of 
Mott insulators in 2 + 1 dimensions P, 13 ■ A closely re- 
lated problem concerns the stability of U{1) spin liquids 
in 2 -f- 1 dimensions [^Q. In either case, models which 
are often considered as toy models in the high-energy 
physics literature are supposed to describe the low-energy 
properties of real systems in condensed matter physics. 
For instance, a model that frequently appears in the 
condensed matter literature is the f2 -I- l)-dimensional 
quantum electrodynamics (QED3) |^. It emerges, 
for instance, as an effective theory for Mott insulators 
0,0, Pi- Let us briefly recall how it arises in this context. 
The Hamiltonian of a SU{N) Heisenberg antiferromag- 
net is written in a slave-fermion representation as H ^ 
-{J/N)J2{t,j) fLfjafjpfifi, where the local constraint 
fiafia = N/2 holds. A Hubbard-Stratonovich transfor- 
mation introduces the auxiliary field Xij — ifiafja) 0- 
The resulting effective theory can be treated as a lattice 
gauge theory, where the gauge field Aij emerges as the 
phase of Xijj i-e-, Xij = Xo^''^'^ , where xo is determined 
from mean-field theory. The (2 + l)-dimensional low- 
energy effective Lagrangian in imaginary time has the 
form ^ r 
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where each ipa is a four-component Dirac spinor and 
Ffj.,^ = dfj,Ai, — duA^ is the usual field strength tensor. 
A rough estimate of the bare gauge coupling is given by 
Cq ~ Xo'^'^j where a is the lattice spacing. 

An anisotropic version of QED3 has also been stud- 
ied in the context of phase fluctuations in rf-wave super- 
conductors 0, Q|. A key feature of the QED3 the- 
ory of Mott insulators is its parity conservation. In 



fact, it is possible to introduce two different QED3s, one 
which conserves parity and one which does not. The lat- 
ter theory involves two-component spinors, and allows 
for a chirally-invariant mass term which is not parity- 
invariant. In such a QED3 theory a Chern-Simons term 
is generated by fluctuations The QED3 the- 

ory relevant to Mott insulators and d-wave supercon- 
ductors involves four-coinponent spinors, and does pos- 
sess chiral symmetry |5|,|6|. In such a model, the chiral 
symmetry can be spontaneously broken through the dy- 
namical generation of a fermion mass. In the context 
of Mott insulators, the chiral symmetry breaking corre- 
sponds to the development of Neel order . Indeed, the 
nonzero condensate {ipip) corresponds to the staggered 
magnetization. Since in condensed matter physics parity- 
conserving QED3 is not just a toy model, and that the 
four-component Dirac spinors represent physical excita- 
tions — the low-energy spinous, we may call this theory 
quantum spinodynamics. 

An important feature of QED3 for Mott insulators is 
that the U{1) gauge group is compact. The compactness 
causes important changes in the physical properties of 
the theory. It allows for quantum excitation of magnetic 
monopoles which play an important role in determining 
the phase structure of the theory. This has been known 
for a long time. In particular, Polyakov 14] has shown 
that compact Maxwell theory in 2 + 1 dimensions con- 
fines permanently electric test charges. The electrostatic 
potential has the form V{R) ~ R, instead of the usual 
two-dimensional Coulomb potential V{R) ~ In i? of the 
non-compact Maxwell theory in 2 -I- 1 dimensions. Since 
V{R) ~ R holds for all values of the gauge coupling, 
the compact (2 + l)-dimensional Maxwell theory does 
not exhibit any phase transition, i.e., the confinement 
is permanent. This theory is equivalent to a Coulomb 
gas of magnetic monopoles in three dimensions and it is 
well known that such a gas does not undergo any phase 
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transition. However, when matter fields are included the 
situation changes, and a deconfinement transition may 
occur. Indeed, matter fields induce shape fluctuations of 
the electric flux tube, leading to a correction term to the 
linearly confining potential. Using a string model for the 
electric flux tube, Liischer 0| found 



V{R) = aR 



(d- 2)7r 
24i? 



(2) 



where a is the string tension. At the deconfining critical 
point, the string tension vanishes and only the Liischer 
term remains at long distances. It represents the "black- 
body" energy of the {d—2) transverse fluctuations of the 
two-dimensional worldsheet of the string . 

Interestingly, by studying QCD near four dimensions, 
Peskin ^3 found that, at the critical point, the in- 
terquark potential does have the 1/i?- behavior for all 
d S (4, 4 -|- e) , and argued that this should also be valid 
outside this small dimension interval. Recently we 
have found that such a behavior is also realized in an 
Abelian U(l)-gauge theory for d G (2,4), provided that 
this is coupled to matter fields. 

In order to better illustrate this mechanism, we shall 
explicitly perform the calculation of the interspinon po- 
tential to one-loop order in arbitrary space-time dimen- 
sion 2 < d < 4, going to the case of interest d = 3 at the 
end 0. The potential is defined by 
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where ep is the bare electric charge, and n(g) is the 
vacuum polarization. At one-loop order, the vaccum 
polarization is given by Ii-{q) = ^A{d)N e'^\q\'^~'^ , with 
A{d) = r(2 - (i/2)r2(d/2)/[(47r)''/2r(d)]. At large dis- 
tances the vacuum polarization gives the more relevant 
contribution to the scaling behavior if 2 < d < 4, and the 
interspinon potential is given by 
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For d — i the above potential becomes simply V{R) = 
—A/{-kNR). Interestingly, by expanding |0J near d — 2, 
we obtain at lowest order 



V{R) 



(d- 2)7r 
8iVi? ' 



(5) 



which has for iV = 3 replica precisely the form of the 
Liischer term, although this theory has no confinement. 
In order to allow for this, we compatify the U(l) gauge 
group which gives rise to magnetic monopoles. In the 
absence of fermions this theory in d = 3 is known to 
be described via a duality transformation by the sine- 
Gordon Lagrangian 
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This is also the field theory of a Coulomb gas of 
monopoles, with zq being the bare fugacity of the gas. 
The RG equations for a Goulomb gas in d = 3 were 
obtained a long time ago by Kosterlitz His re- 

sults can be used here to obtain the RG equation for 
the gauge coupling in the absence of fermionic matter. 



By introducing the dimensionless couplings / ; 



and y = z{l)/{eQ)^, where I — In(eQr) is a logarithmic 
length scale, we obtain 
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The above equations imply that there is no fixed point 
for the gauge coupling. Therefore, the compact three- 
dimensional Maxwell theory does not undergo any phase 
transition. The photon mass = Stt^z/c^ is always 
non-zero and the theory confines permanently the electric 
charges. It can be seen from Eq. that a kind of anti- 
screening happens in this theory, which is responsible for 
confinement. Indeed, we can rewrite Eq. Q a,s df /dl — 
(1 — 7a)/, with 7A = — 47r^y^//. The negative sign of 
jA is actually a remarkable example of the intimate link 
between asymptotic freedom and confinement . 

Next we obtain the modification to Eqs. iQ and 
(O due to the coupling with the matter fields. In or- 
der to derive the RG equations including matter we 
have employed a formalism similar to the one devel- 
oped by Young in the case of the two-dimensional 
Coulomb gas. This formalism is based on a mean- field 
self-consistent approximation and applies very well to the 
d > 2 case, since d = 2 is the upper critical dimension 
for the Coulomb gas. The needed modification comes 
from the extra renormalization of the gauge coupling due 
to the vacuum polarization. This leads to an effective 
charge e^(/) = e(e'/e§)ZA(Oe'eQ, where e{r) is the scale- 
dependent "dielectric" constant of the Coulomb gas of 
magnetic monopoles, and Za{1) is the gauge field wave 
function renormalization. In terms of dimensionless cou- 
plings this leads to 
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where 7^ = — dlnZ^/dL However, the expression of 
in terms of bare variables is different from before, be- 
ing given by y^ = (327^V3)zgZA(0e^'""('V(eo)^ where 
u{l) = [/(e'/cg) is a self-consistent magnetic monopole 
potential satisfying du/dl = i^'^ / f [l^. Therefore, the 
coupling to matter modify also Eq. (|SJ) to 



rfy ^ /„ _ TT^ _ 7A 
dl \ f 2 



y- 



(10) 



Note the crucial difference between the analysis made 
here and the one of Refs. and (2^. There it was 
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assumed that the underlying non- compact theory is crit- 
ical, and monopoles were introduced only at that point. 
This corresponds to take the RG function o-t the fixed 
point of the non-compact theory, i.e., 7^ = 77^ = 1 [2^ 
for d — 3. In this way, Eqs. © and l(Tn|l become similar 
to the RG equations of a Kosterhtz-Thouless phase tran- 
sition j27l | , except that the present dimensionahty is three 
instead of two Jlliil. Our Eqs. © and ^ have the 
advantage of being vahd at all length scales. Eqs. © and 

rare similar to the ones in the work of Hermclc et al. 
There is, however, an important difference: Eq. H1U|I 
contains the correction proportional to which is absent 
in Ref . 3 ■ This will allow us to strenghten considerably 
the results obtained by these authors. It is important to 
emphasize that the additional term in Eq. 1)1 0|l cannot 
be neglected even if a large N limit is assumed. Indeed, 
since the large N limit is taken for Nbq fixed, it follows 
that J A ^ as it should be, since it gives the anoma- 

lous dimension RG function of the non-compact theory. 
It is of the same order as the first term between parenthe- 
ses in Eq. 1)10(1 . which corresponds to the dimensionality 
of the space-time. Thus, this problem has no obvious 
control parameter, and should be seen as a matter field 
fluctuation-corrected Debye-Hiickel theory. In the case of 
compact Maxwell theory, the Debye-Hiickel theory corre- 
sponds to a non-dilute gas of monopoles and its validity is 
determined the parameter nX^ , where n is the monopole 
density and Xjy = \J e^ jiA.-K^ri) is the Debye Length. For 
the compact Maxwell theory we have that nA^j ^ 1 
and the Debye-Hiickel theory is a very good approxima- 
tion. Including matter fields makes the monopole gas 
dilute and nA|j is no longer large. The Debye-Hiickel pa- 
rameter can be written as riA|) — '\/2e^/(87r^M), where 
M = 2-K\f2zle is the photon mass. Thus, in the pres- 
ence of matter a perturbation theory around the compact 
Maxwell theory can be performed where e^/M is a small 
parameter. We will see below that the fixed points at 
nonzero fugacity give indeed a small value of nA|). 

By considering the one- loop result 7^ — N f we 
find besides the fixed points /* = and = of 

the non-compact theory, the following non-trivial fixed 
points governing the phase structure of compact QED3: 



/± = 1 (e ± ^36 - Nti^ 



y± 



1 / 60 - Nn^ ± 10\/36 - Nn^ 

TT \ 2N 
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The above fixed points exist only for N < Nc = 36 / tt^ ^ 
1.161. For N > Nc only the non-compact fixed point 
exists. For = 1 the fixed points {f±,y±) give for 
the Debye-Hiickel parameter the values (nA|j)+ « 0.3 
and (nA|j)_ « 0.15, respectively. In Fig. 1 we show 
a schematic flow diagram for the case = 1. The 
dashed line in the flow diagram passes through the unsta- 
ble fixed point having coordinates (/_ , y- ) . This line sep- 



FIG. 1: Schematic flow diagram for the case A'^ = 1 



arates two different critical regimes. Note that 7^1 is A^- 
dependent at the fixed points {f±,y±), with its two pos- 
sible values given by 7^ = A^/±/8 = (6± V36^^^77^)/2. 
This result implies that there is no N > 1 for which 
7^ = 1 in compact QED3. This rules out a KT-like 
transition in compact QED3 for physical values of N. 

The flow diagram in Fig. 1 indicates two distinct phys- 
ical regimes governed by stable fixed points separated by 
the dashed line in the figure. Depending on the initial 
conditions on the physical parameters, the system will 
choose to flow either to the non-compact fixed point be- 
low the dashed line, or to the compact one above the 
dashed line. The interesting physical regime for us is 
governed by the fixed points at nonzero fugacity. It is 
clear that the fixed points {f±,y±) are associated with 
confined phases, since there both the photon mass M 
and string tension a = 2e^M/7r^ are nonzero. The string 
tension approaches a universal value as TV approaches Nc 
from the left, i.e., lim^y^^- a/ef^ — 8(7r/3)'^/^. Since for 
N > Nc the string tension vanishes, it follows that there 
is a universal jump at Nc. Thus, in the present con- 
text the string stiffness behaves similarly to the super- 
fluid stiffness in two-dimensional superfluids '29*1 , though 
here there is no KT transition. The vanishing of the 
string tension above Nc is a clear signature for spinon 
deconfinement for N = 2. 

Below Nc the interspinon potential has the form 
V{R) = o-i? - a/R + 0(l/i?2), where a is the univer- 
sal coefficient of the Liischer term for the string fluctu- 
ation in compact QED3. The coefficient a is deflned 
hy a = /c/27r, where fc is any of the three charged 
fixed points in Fig. 1. For the stable confining regime 
governed by the fixed point (/-|_, ?/-(_) we obtain that 
a 2(6 + V36- 7V7r3)/7riV. 

It is perfectly plausible to argue that in the confined 
phase of compact QED3 the chiral symmetry is broken, 
just as in the QCD case ■ Chiral symmetry break- 
ing is believed to occur in QED3 for N < Nch, where 
typically A'ch 3. Indeed, an early estimate based 
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on the analysis of the Schwinger-Dyson equation gives 
iVch = « 3.2 0, which was roughly confirmed 

by a Monte Carlo simulation giving A^ch = 3-5 ± 0.5 [sij . 
However, the true value of N^h is still far from being con- 
sensual. For instance, recent Monte Carlo simulations 
do not find a decisive indication that chiral symmetry is 
broken for N > 2 [s^ and an elaborate aiialysis of the 
Schwinger-Dyson equations gives A^ch « 4 PJ. Our re- 
sults indicate that for Nc < N < Nch the spinous are 
deconfined but chiral symmetry is broken. However, it 
is not excluded that Nc = Nch- Recently, a conjectured 
inequality was used to suggest that iVch — 3/2 in the 
non-compact case j^. It is remarkable that our critical 
value Nc is so close to the latter estimate. If Nch > 2, we 
would obtain that for the physical case N — 2 antifer- 
romagnetism is present 9] , while the spinous are decon- 
fined. In such a situation doping will eventually destroy 
the magnetic order and, since the spinous are deconfined, 
a genuine spin liquid will develop. Our results confirm 
the analysis of Ref. 0], whose discussion was made in 
the large N limit. 

We are indebted to Zlatko Tesanovic, who pointed out 
an important mistake in a previous version of the pa- 
per. This work received partial support of the european 
network COSLAB. 
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